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The WT1 gene encodes a zinc finger transcription factor important for normal kidney development. WT1 is a sup- 
pressor for Wilms tumour development and an oncogene for diverse malignant tumours. We recently established 
cell linesfrom primary Wilms tumours with different WT1 mutations. To investigate thef unction of mutant WT1 pro- 
teins, we performed WT1 knockdown experiments in cell lines with a f rameshift/extension (p.V432fsX87 = Wilms3) 
and a stop mutation (p.P362X = Wilms2) of WT1, followed by genome-wide gene expression analysis. We also 
expressed wild-typeand mutant WT1 proteins in human mesenchymal stem cellsand established gene expression 
profiles. A detailed analysis of gene expression data enabled us to classify the WT1 mutations as gain-of-function 
mutations. The mutant WT1 Wllms2 and WT1 Wllms3 proteins acquired an ability to modulate the expression of a highly 
significant number of genes from the G2/M phase of the cell cycle, and WT1 knockdown experiments showed that 
they are required for Wilms tumour cell proliferation. p53 negatively regulates the activity of a large number of these 
genes that are also part of a core proliferation cluster in diverse human cancers. Our data strongly suggest that 
mutant WT1 proteins facilitate expression of these cell cycle genes by antagonizing transcriptional repression 
mediated by p53. We show that mutant WT1 can physically interact with p53. Together the findings show for the 
first time that mutant WT1 proteins have a gain-of-function and act as oncogenes for Wilms tumour development 
by regulating Wilms tumour cell proliferation. 



INTRODUCTION 

Wilms tumour is a paediatric kidney cancer affecting 1/10 000 
children a year. The first protein to be associated with WT devel- 
opment is encoded by the gene WT1 located on chromosome 
1 lpl3 (1,2). WT1 is mutated in ~ 15-20% of all WT and is an 
important factor for normal kidney development (3). The gene 
encodes a protein of 52-54 kDa with exons 7 to 10 encoding 
four C 2 -H 2 zinc fingers (ZFs) of the Kriippel type that bind 
DNA and RNA. The first exons encode a proline -glutamine 



(Pro/Gln)-rich domain that contains a putative RNA recognition 
motif and is involved in transcriptional repression and activa- 
tion, dimerization and nuclear localization (4-7). 

Alternative splicing results in four major isoforms, the first 
leading to inclusion/exclusion of exon 5 and the second to inclu- 
sion/exclusion of three amino acids, lysine, threonine and serine 
(KTS) after exon 9. It was first shown that WT1 lacking KTS 
binds to a GC-rich EGR1 consensus sequence, as well as to 
an unrelated TCC repeat motif (8,9). The inclusion of KTS 
between ZF3 and 4 significantly reduces the DNA-binding 
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affinity of WT1 and the +KTS isoform binds to other DNA 
targets (10). There is also evidence that both WT1 isoforms + 
and — KTS are involved in post-transcriptional processes (11). 
The +KTS isoform co-localizes and co-immunoprecipitates 
with splice factors, and WT 1 can modify splicing by interacting 
with the splice factor U2AF65 (12,13). Using the RNA selection 
method SELEX and WT 1 ZF constructs, three RNA aptamers 
that are recognized by WT1 were identified (14). Three of four 
ZFs were necessary, and deletion of ZF1 resulted in reduced 
and insertion of KTS abolished binding for the RNA targets 
(14,15). Using these RNA aptamers, Weiss and Romaniuk 
showed that ZF2 and 3 are necessary for RNA binding (16). 

WT 1 was found in poly(A) + nuclear RNP from foetal kidneys 
(17) and in mRNP particles in K562 cells, further pointing to 
a role in post-transcriptional regulation. There is also strong 
evidence that WT1 binds to mRNA in vivo with an important 
role of ZF1 in RNA binding (17). a-Actininl (ACTN1) was 
identified as the first cellular RNA target for WT1 (18). In add- 
ition, WT1 is also found in the cytoplasm, specifically in 
cancer cells, suggesting that WT1 could be involved in mRNA 
export, translation or turnover. 

WT 1 can form homo- and hetero-dimers with various other pro- 
teins and is found in large complexes in Wilms tumours (19,20). 
WT 1 can activate or repress genes depending on the architecture 
of the promoter; the transcriptional function of WT1 can be 
modulated through the physical interaction with other proteins 
(20). The N-terminal transcriptional regulation domain of WT1 
was found to bind BASP1, resulting in transcriptional co- 
suppressor activity (21). WT1 and p53 co-immunoprecipitated 
in adenovirus ElA-immortalized baby rat kidney cells, and this 
association resulted in stabilization and increased levels of 
p53 with enhanced binding to and transactivation of a promoter 
containing the p53 RGC target sequence. In contrast, repression 
of a TATA-containing promoter by p53 was inhibited by WT1 
(22). The WT1 domain, necessary for binding to p53, was 
mapped to ZF1 and 2. 

Most mutations identified in Wilms tumours lead to a trun- 
cated WT1 protein that has lost all or parts of the ZF domain, 
suggesting that these proteins are impaired in binding to the 
EGR consensus or the alternative TCC WT 1 -binding motif. In 
most cases, however, the N-terminal Pro/Gln-rich repression 
and activation domain is intact and it is possible that the 
mutant proteins can still interact with other proteins, thereby 
modulating their activity. Such a scenario where a truncated 
mutant WT 1 protein inhibits the transcriptional activity of wild- 
type WT1 has been described (23). In addition, truncated mutant 
WT 1 proteins alter the subcellular localization of the wild-type 
protein (24). In some cases, the mutations lead to truncation in 
or after ZF2, 3 or 4, and these proteins may also have retained 
their RNA-binding capacity or they may bind altered DNA 
motifs. Recently, a nuclear localization signal in the third ZF 
was identified that is required for nuclear import (25). We have 
shown previously that mutant WT 1 proteins are predominantly 
localized to the cytoplasm in Wilms tumours (26) and in 
Wilms tumour cells as well as in NIH3T3 cells after transfection 
with mutant WT1 (27). 

We have previously described a method for the successful es- 
tablishment of Wilms tumour cell lines from Wilms tumours 
with WT1 mutations (27). All cell lines carry a homozygous 
WT1 mutation owing to loss of heterozygosity of 1 lp markers. 



Only one cell line from a WAGR patient has a hemizygous 
WT1 mutation on the remaining allele (Wilms4). These cell 
lines can be grown for ~ 20 passages but do not have an unlimited 
life span. With this unique Wilms tumour cell culture model 
system, where both alleles of WT1 are mutant and no wild-type 
WT1 allele is present, we can now begin to study for the first 
time the function of the mutant WT 1 proteins in a homologous 
system (27). We have previously shown that the Wilms2 cell 
line has a WT1 stop mutation in exon 8 leading to a truncation 
in ZF2 (p.R362X= WTl wilms2 ) and a p.S45Y mutation in 
CTNNB1. The WT1 frameshift mutation in exon 10 of the 
Wilms3 tumour cell line leads to an elongation of the WT1 
protein by 68 amino acids (p.V432fsX87 = WTl Wilms3 ); this 
cell line is wild type for CTNNB1. In Wilms3 cells, the mutant 
RNA is expressed, the WT1 protein is present in the cytoplasm 
of the cells, and all isoforms of the WTl Wllms3 mutant proteins 
are localized in the cytoplasm of transfected NIH3T3 cells 
(27). Both WT1 mutant proteins retain the Pro/Gln-rich 
N-terminus and an intact ZF1. In the WTl wilms3 mutant 
protein, only ZF4 is disrupted, but both + and — KTS isoforms 
can be synthesized, whereas in Wilms2, ZF2 is truncated and 
ZF3 and 4 are missing and the stop codon occurs before the 
KTS alternative splice site. The aim of this study was to function- 
ally characterize these two naturally occurring mutant WT 1 pro- 
teins and to establish whether they act via gain or loss of function. 
Here, we report that two cytoplasmically expressed mutant WT 1 
proteins regulate a number of genes/proteins as revealed by 
knockdown of mutant WT1 in the Wilms tumour cell lines. 
In this work, we show that the mutant proteins retain their 
ability to interact with p53 and to bind to RNA with a reduced 
association constant. Loss of WT1 by knockdown in these cells 
results in reduced proliferation and reduced expression of 
genes from the G2/M phase of the cell cycle. Expression of the 
mutant WTl Wllms3 protein in mesenchymal stem cells (MSCs) 
results in up-regulation of the same cell cycle genes, and 
these genes are not regulated by wild-type WT1, confirming a 
gain-of-function for the mutant protein. 



RESULTS 

Intracellular distribution of mutant WTl Wilms2 protein 

Wilms tumour protein WT 1 is localized to the nucleus in embry- 
onic and adult tissues (28). We have reported recently that 
contrary to wild-type WT1, mutant WTl-GFP-fusion proteins 
from Wilms3 tumour cells are located predominantly in the 
cytoplasm of transfected NIH3T3 cells (27). The mutation in 
^ywiimss leads tQ an extended WT1 p ro tein with additional 

68 amino acids at the C-terminus. The positions of mutations in 



WTl wiims3 and WTl Wllms2 are shown in Figure 1. WTl Wllms2 is 
a truncated protein, and the stop codon occurs before the KTS 
alternative splice site. Therefore, only the isoforms + and - 
exon 5 exist. In both Wilms tumour cell lines, the WT1 mutation 
is homozygous, i.e. both alleles carry the same mutation and no 
normal WT1 protein is present in the cells. Here we investigated 
the intracellular distribution of these proteins in NIH3T3 cells. 



WT1 Wilms2+exon5_ WTJ 



Wilms2— exon5 



GFP-fusion proteins 



were expressed in NIH3T3 cells, and a predominant cytoplasmic 
localization was observed by fluorescence microscopy (Supple- 
mentary, Material Fig. SI). In contrast, the wild-type WTl~ /_ 
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and _/+ isoforms are exclusively located in the nucleus (Supple- 
mentary Material, Fig. SI). Taken together, these results show 
an important difference between wild- type and mutant WT1 
proteins. 



Interaction of mutant WT1 proteins with p53 

It has been reported that wild-type WT 1 protein has the ability to 
interact with the tumour suppressor p53 (22). Here we wished to 
examine whether mutant WTl Wilms2 and WTl Wilms3 proteins 
have retained this ability. To investigate whether mutant WT1 
proteins interact with p53, HeLa cells that express p53 but no 
WT 1 protein were transfected with plasmids driving expression 
of GFP, or GFP linked to wild-type WT1, WTl wrfms2 and 
WTl Wllms3 . Preparation of whole cell extracts, followed by 
immunoprecipitation with antibodies specific for p53 or a 
control antibody and immunoblotting with a WT1 antibody spe- 
cific for the N-terminus, demonstrated co-immunoprecipitation 
of these WT1 proteins with p53. As reported before, wild-type 
WT1 co-immunoprecipitated with p53 (Fig. 2). Moreover, 
both WTl Wilms2 and WTl Wilms3 retained the ability to interact 



with p53 (Fig. 2). It is thus possible that mutant WT1 proteins 
modulate p53 activities in Wilms tumour cells. 
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Figure 1. Schematic drawings of wild-type and mutant WT1 proteins. Wild-type 
WT 1 protein (top). Red asterisks indicate the positions of WTl mutations that we 
identified previously in Wilms2 and Wilms3 tumour cells. Truncated WT 1 v 



protein (middle). Extended WT1 



(bottom). The dotted red flag at the 



Equilibrium binding of RNA aptamers to the ZF 
domains of WT1 and Wilms3 

As the mutant WT 1 Wllms3 protein harbours an intact ZF 1 , 2 and 3 , 
we wished to study whether the ZF domain from this protein 
retains its capacity to bind to RNA. The RNA-binding activity 
of the WT1 ZF domain has been characterized previously 
using the RNA aptamers RNA22 and RNA38 (29). Here we 
have determined the affinities of the Wilms3 ZF domains for 
these two different aptamers using an equilibrium binding 
assay (Fig. 3). Although the mutation in Wilms3 disrupts 
proper zinc coordination to finger 4 of the domain, RNA-binding 
activity is retained but reduced by 4- to 5-fold compared with the 
wild-type WT1 ZF domain (Fig. 3). 



Mutant WT1 proteins are linked to gene regulation 
in Wilms tumour cells 

To gain more insie ht into the function of mutant WT1 wilms3 and 
WT1 wums2 protemS) we p er f ormec i s iRNA (siWTl) knockdown 

analyses. First, we tested different ffTV-specific siRNA oligonu- 
cleotides and found that a combination of two siRNAs directed 
against sequences in exon 6 and the 3' UTR of WT1 achieved a 
very efficient knockdown 48 h after transfection of Wilms3 
cells (Fig. 4A). As control a non-silencing siRNA was used. 
Gene expression profiles were generated using whole-genome 
microarrays. To identify WT1 responsive genes, the mean of 
two biological replicates was determined and the corresponding 
genes were selected (P < 0.05). A 1.5-fold change (fc) relative 
to the control and a minimal expression of 200 fluorescence in- 
tensity was used as cut-off, and we detected 388 genes with 
increased and 381 genes with decreased expression. The effi- 
ciency of WTl mims3 knockdown was examined by 
Q-RT-PCR, and a reduction of WTl Wilms3 expression by 82% 
was measured (Fig. 4A). Microarray data were validated for 
selected genes by Q-RT-PCR (Fig. 4B). It is evident that 
siWTl knockdown caused decreased expression of GDF5 and 
increased expression of MESCD2, CASP7, PODXL1 in 
Wilms3 tumour cells. We next examined gene expression pro- 
files after siWTl knockdown in Wilms2 cells, and the efficiency 
of WTl" 



Wilms2 



C-terminus of WTl lms " indicates the position of the changed 87 amino acids. 



knockdown was quantified by Q-RT-PCR 
(Fig. 4C). We identified 361 genes with increased and 422 
genes with decreased expression levels. Microarray data 
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Figure 2. Interaction of p53 and WT 1 proteins. Detection of WT 1 -p53 protein interactions by immunoprecipitation and western blot analysis. HeLa cells were trans- 
fected with expression plasmids encoding the WTl-GFP-fusion proteins, wild-type WTl, mutant WTl Wllms3 and mutant WTl Wllms2 . As control, cells were trans- 
fected with a GFP-expression plasmid alone. Lysates from transfected cells were incubated with control antibodies (lanes, aCtrl) or p53-specific antibodies 
(lanes, aP53). Immunoprecipitates were resolved by SDS-PAGE and immunoblotted using antibodies specific for the N-terminus of the WTl protein 
(aF6WT 1 ). As input controls, 5% of the respective lysates from transfected HeLa cells were loaded on the gel (lanes, In). Asterisk indicates a non-specific band. Wild- 
type and mutant WTl-GFP-fusion proteins in the input and immunoprecipitated fractions are marked by black boxes. 
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Figure 3. Equilibrium binding curves for the binding of the ZF domains of WT 1 
and Wilms3 to RNA. The affinities of RNA22 and RNA3 8 were measured using a 
nitrocellulose filter binding assay. Each line represents the best fit for the forma- 
tion of a simple bimolecular complex (R values of 0.98-0.99). Each data point is 
the mean of two independent determinations, with the standard deviation indi- 
cated with error bars. The table shows the apparent association constant for wild- 
type WT1 and mutant WTl wilms3 . 



showed a 2.7-fold up-regulation of PODXL1, which was vali- 
dated by Q-RT-PCR (Fig. 4C). Next we compared the si WT1 
gene expression data from Wilms2 and Wilms3 tumour cells 
and identified a common set of 280 down- and 243 up-regulated 
genes (Supplementary Material, Fig. S2). These results demon- 
strate that the mutant WT1 proteins are engaged in regulating 
gene expression in Wilms tumour cells. Moreover, the data 
show that truncated WTl Wllms2 and extended WTl Wllms3 pro- 
teins (Fig. 1) regulate similar gene expression programmes. 

As this result is of great importance for the field of WT1 
research, we wished to confirm our findings by another 
method. We thus performed lenti viral-mediated shWTl knock- 
down experiments using Wilms3 tumour cells. The shWTl 
targets exon 7 of the WT1 gene, and the most efficient knock- 
down of WTl mlms3 expression (73%) occurred 48 h after infec- 
tion (Fig. 4D). Gene expression analyses were performed under 
these conditions, and a scrambled non-silencing shRNA was 
used as control. A 1 .5 fc relative to the control and a minimal ex- 
pression of 200 fluorescence intensity was used as cut-off. The 
microarrays identified 1012 genes with increased and 1076 
genes with decreased expression levels. Increased expression 
of PODXL and JARID2 was validated by Q-RT-PCR 
(Fig. 4E). Taken together our WT1 knockdown experiments 
show for the first time that mutant WT 1 proteins do have a bio- 
logical function and this is linked to gene regulation in Wilms 
tumour cells. The top 100 genes with decreased and increased 
expression are listed in Supplementary Material, Tables SI 



and S2, showing the fold changes of expression relative to the 
scrambled control. 

Mutant WT1 proteins promote expression of genes 
from the G2/M phase of the cell cycle 

To gain more insight into the biological function of the mutant 
WTl Wllms3 protein, we analysed the gene expression data from 
siWTl knockdown experiments using MetaCore, an integrated 
software suite for functional analysis of biological experimental 
data. To identify biological processes associated with gene 
expression data, we used an enrichment analysis of MetaCore 
pathway maps that represent a set of signalling and metabolic 
maps covering cellular processes in a comprehensive way. 
These maps also show functions and interactions of individual 
genes in a given pathway. We first examined the subset of down- 
regulated genes obtained from siWTl experiments in Wilms3 
cells using MetaCore for enrichment of specific pathways and 
identified highly significant enrichment of several maps from 
the cell cycle. A list of the ten most significant pathway maps 
and their P-values is summarized in Figure 5A (blue bars). The 
genes from the cell cycle pathways are listed in Supplementary 
Material, Table S3, left column. It is evident that these pathways 
are related to the G2/M phase of cell cycle. The most significant 
pathway maps ' Cell cycle_Role of APC in cell cycle regulation' , 
'Cell cycle_Spindle assembly and chromosome separation' and 
'Cell cycle_The metaphase checkpoint' operate at the end of G2/ 
entry into mitosis. 

Next, we analysed the subset of genes that were down-regulated 
after infection of Wilms3 cells with shWTl lentiviruses. The Venn 
diagram in Supplementary Material, Fig. S3 Panel A, shows that 
70 genes are regulated in common in shWTl and siWTl knock- 
down experiments. Enrichment analysis of MetaCore pathway 
maps revealed 'Cell cycle_The metaphase checkpoint', 'Cell 
cycle_Chromosome condensation in prometaphase' and 'Cell 
cycle_initiation of mitosis as highly significant' . A list of the ten 
top-ranking pathway maps is shown in Figure 5A, orange bars. 
The most significant pathway map 'Cell cycle_The metaphase 
checkpoint' is visualized in Figure 5B, and the down-regulated 
genes are indicated by blue thermometers. The genes correspond- 
ing to the cell cycle pathways are listed in Supplementary Material 
Table S3, right column. Although only 70 identical genes are regu- 
lated by si and shRNA, Supplementary Material, Table S3 shows 
that additional but different genes from the same pathways are 
regulated in either approach. These data show for the first time 
that mutant WT 1 from Wilms3 tumour cells promote the expres- 
sion of a significant number of genes from the G2/M phase of the 
cell cycle. Pathway maps, however, show only selected genes and 
do not represent all genes of a given biological process. To get a 
comprehensive overview, we performed a GO process analysis 
using MetaCore and the subset of down-regulated genes. We iden- 
tified the most significant biological processes associated with the 
set of down-regulated genes (Supplementary Material, Fig. S4). 
These included 'cell cycle' (P-value 1.113E— 15), 'mitosis' 
(P-value 1.5 HE- 14) and 'cell division' (P-value 5.459E- 14). 

We next performed a MetaCore pathway enrichment analysis 
of the subset of genes whose expression was down-regulated 
after siRNA-mediated knockdown of mutant WTl Wilms2 . The 
enrichment analysis of MetaCore pathway maps revealed 
that 'Cell cycle_Role of APC in cell cycle regulation', 'DNA 
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Figure 4. Validation of microarray data from si WT1 and sh WT1 knockdown experiments by Q-RT-PCR. All Q-RT-PCR data are displayed relative to a non-silencing 
si- or shRNA control that is set at 1 . Data represent the mean of triplicates, and the standard deviations are indicated for each experiment. PCR reactions were done after 
transfection of siRNA, and the time points are indicated for each Q-RT-PCR. (A) WT1 w,ll " s3 expression after treatment with si W 77-RNA. (B) Validation of GDF5, 



CASP7. MESCD2 and PODXL expression after siRNA-mediated WT1 knockdown in Wilms3 cells. (C) WT1" 



and PODXL expression after treatment with 



siffTZ-RNA. (D) WTl w ' l "' i3 expression 48 h after infection with shWTl lentivimses (left panel). Expression data are displayed relative to a non-silencing shRNA- 
lentivirus control that is set at 1 . Validation of JARLD2 and PODXL1 expression in Wilms3 cells infected with shWTl lentiviruses (right panel). Q-RT-PCR was per- 
formed 48 h after infection with sbWTl lentiviruses. 



damage_ATM/ATR regulation of Gl/S checkpoint' and 'Cell 
cycle_the metaphase checkpoint', are the three most significant 
pathway maps. A list of the ten top-ranking pathway maps and 
their P-values for enrichment is shown in Figure 5C, and the 
corresponding genes are shown in Supplementary Material, 
Table S4. We can conclude here that the functions of truncated 

WT1 Wilm S 2 an( j extended WT iWilms3 proteins are very s i m il ar 

in the context of gene regulation. 

WTl Wllms3 inhibits the expression of genes from multiple 
signalling pathways 

Next, we wished to examine the biological function of genes 
that are repressed by mutant WTl Wllms3 . We therefore used the 
subsets of genes that were up-regulated in siWTl and shWTl 
knockdown experiments and performed enrichment analyses of 
pathway maps using MetaCore. We noticed that the P-values for 
pathway enrichment were lower in siWTl than in shWTl 



experiments and selected the shWTl data from Wilms3 cells for 
our analysis. Enrichment analysis by pathway maps identified 
'Transcription_Androgen Receptor nuclear signalling' and Cyto- 
skeleton remodelling_TGF, WNT and cytoskeleton remodelling' 
as most significant. A list of the top-ranking pathway maps is 
presented in Supplementary Material, Fig. S5 and Table S5. This 
analysis shows an overrepresentation of enriched pathway maps 
related to signalling, and we decided to investigate this issue in a 
more comprehensive fashion. We therefore inspected all 
pathway maps from the MetaCore database related to signalling 
and identified a large number of genes that are repressed 
by mutant WT1. These include TGFBRII, IGF1R, AXL, MET, 
WNT5 and FZD6 and the downstream signalling components 
DSH, CTNNB1, AKT, DUSP3, DUSP16, ERK2, FAK1, MAPK, 
MLCK, MYLK1, STAT2/3, FOS, SOS and others. The proteins 
encoded by DUSP3 and DUSP16 are involved in the negative 
regulation of Erk signalling and MAPK pathways, respectively 
(Supplementary Material, Fig. S6). Moreover, WTl Wilms3 seems 
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Figure 5. Mutant WT1 proteins control expression of cell cycle genes from the 
G2/M phase in Wilms2 and Wilms3 tumour cells. (A) Enrichment of MetaCore 
pathway maps in gene expression profiles from Wilms tumour cells after WTI 
knockdown. For enrichment analysis, the subsets of down-regulated genes 
were used. The height of the blue and orange bars indicates the P- values of enrich- 
ment for individual pathway maps. Orange bars: WTI knockdown in Wilms3 
cells after infection with shWTl lentivirases. Blue bars: WTI knockdown in 
Wilms3 cells after transfection of siRNA. (B) MetaCore pathway map 'Cell 
cycle_The metaphase checkpoint'. Blue thermometers indicate down-regulated 
genes. The numbers on top of each thermometer refers to individual WTI knock- 
down experiments. 1, shWTl Wilms3 cells; 2, siWTl Wilms3 cells. The proper- 
ties of genes from this pathway are indicated by coloured symbols. These 



to repress expression of CDKN1B (p27) and CDKN2B (p 1 5) genes, 
negative regulators of progression from the G 1 to S phase transition 
of the cell cycle. A comparison of the siWTl with shWTl knock- 
down experiments is shown in Supplementary Material, Fig. S3, 
Panel B . The top-ranking pathway ' signal transduction Erk interac - 
tions: inhibition of Erk' is identified by the separate analysis of the 
70 identical genes that are regulated in both approaches. 

Most Wilms tumours with WTI mutations also have mutations 
in CTNNB1 and in consequence an activated Wnt signalling 
pathway. However, Wilms3 cells do not carry a CTNNB1 mutation, 
and the Wnt/(3-catenin pathway is not active as measured by TOP/ 
FOP reporter plasmids (27). To check whether WTl^ 1 " 183 is 
involved in regulating expression of genes from Wnt signalling 
pathways, we used the subset of genes with increased expression 
in si WTI and sh WTI experiments and examined the Wnt signalling 
pathway maps of the MetaCore database. Our analysis revealed 
that WTl Wims3 negatively regulates the expression of several 
genes from the Wnt signalling pathways including WNT5, 
DKK1, FZD6, FZD8, MESDC2, DSH (DVL3), CTNNB1, as well 
as the downstream target VIM (Supplementary Material, 
Fig. S7). It is known that the regulation of this pathway is highly 
complex, with several negative regulators being up-regulated in a 
cell with active Wnt signalling. The up-regulation of ligands, 
receptors and downstream signalling components upon WTI 
removal suggests that the mutant WTI protein may be a negative 
regulator of the Wnt signalling pathway. A detailed analysis of 
this issue will be published elsewhere. 



Mutant WTI proteins promote proliferation of Wilms 
tumour cells 

Our data show that mutant WT 1 proteins are involved in regulat- 
ing the activities of several genes from the G2/M phase of the cell 
cycle. We therefore investigated whether mutant WTI proteins 
also affect the proliferation of Wilms tumour cells. To this end, 
we determined the proliferation rates of Wilms2 and Wilms3 
cells after WTI knockdown using shWTl lentivirases. Figure 6 
shows that the sliWTl knockdown significantly reduces the 
proliferation of both Wilms2 (P = 0.037) and Wilms3 cells 
(P = 0.014) in comparison with control cells transduced with a 
non-silencing shRNA. A highly significant inhibition of cell pro- 
liferation is also visible in tissue culture dishes 72 and 96 h after 
infection of Wilms3 cells with shWTl lentivirases in comparison 
with the non-silencing control (Supplementary Material, Fig. S8, 
A). Next, we wished to verify these results by a WST proliferation 
assay (Supplementary Material, Fig. S8, B). In this assay, 
cell survival is analysed and we found significantly reduced 
cells numbers in shWTl expressing cells. This suggests 
that some cells do not survive without the presence of 
WT 1 Wllms3 . We conclude from these results that the proliferation 
of Wilms2 and Wilms3 tumour cells is dependent on the activity of 
mutant WTl wilms2 and WTl wilms3 proteins. We report here that 

symbols are explained by a MetaCore quick reference guide (Supplementary Ma- 
terial, Table S15). Green arrows indicate positive interactions (activation) and 
red arrows, negative interactions (inhibition). (C) Enrichment of MetaCore 
pathway maps in gene expression profiles from Wilms2 cells after siWTl knock- 
down. For enrichment analysis, the subsets of down-regulated genes were used. 
The height of the orange bars indicates the P-values for enrichment of individual 
pathway maps. 
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Figure 6. Mutant WT 1 proteins promote proliferation of Wilms tumour cells. Wilms2 and Wilms3 cell number after sh WT1 and scrambled RNA (shC) knockdown at 
different time points, performed in two independent analyses (a and b), showing a significantly lower median cell count after shWTl knockdown compared with 
scrambled RNA, Wilms2: P = 0.037 (Wilcoxon test) and Wilms3 P = 0.014 (Wilcoxon test). 



mutant WT1 proteins facilitate the expression of G2/M phase 
genes, and we asked whether the distribution of cell cycle 
phases would be altered in Wilms3 tumour cells after shWTl 
knockdown. To this end, we used the Cycletest™ Plus kit in 
initial experiments and analysed Wilms3 cells 48 h after infection 
with shlVTl lentiviruses using a fluorescence-activated cell sorter 
(Supplementary Material, Fig. S9). At this time point, we did 
observe a slight increase of Wilms3 cells in G2/M phase in com- 
parison with the controls; however, apoptosis was not observed in 
these studies. Future experiments will determine the kinetics of 
siWTl and shWTl -mediated knockdowns in Wilms3 cells and 
correlate these data with an analysis of cell cycle phases in 
Wilms3 tumour cells. 



WTl wllmsJ protein modulates the activities of selected 
transcription factors including p53 

Our data show that mutant WT 1 proteins are linked to gene regu- 
lation. As both Wilms2 and Wilms3 mutant proteins lack a func- 
tional DNA binding domain, it was of interest to gain more 
insight into the underlying molecular mechanisms. To investi- 
gate this issue, we used the MetaCore program 'What are the 
key transcription factors and target genes in my data?' This 
program uses a database containing all published information 
regarding transcription factors and their target genes, and the 
results show which transcription factors play a significant role 
in individual microarray experiments. First we uploaded the 
down-regulated genes from shWTl Wilms3 knockdown experi- 
ments as the input list. Using these data, the transcription regula- 
tion workflow generates sub-networks centred on transcription 
factors. Sub-networks are ranked by a P-value and interpreted 
in terms of Gene Ontology. The resulted sub-networks of the 
whole network are displayed in Supplementary Material, 
Table S6, a network list table, in the order reflecting their statis- 
tical significance in terms of P-values. The network list shows 
the top 21 most significant sub-networks centred around 
CPvEBl, c-MYC, SP1, YY1, OCT 3-4, GCR-alpha, ESR1, 
p53, E2F1, Androgen Receptor, HIF1A and others. The 



column 'P-value for enriched transcription factor' indicates 
the level of significance. The column 'enriched GO processes' 
lists specific GO processes that were identified in each transcrip- 
tion factor sub-network. For example, the sub-networks centred 
around CREB 1 , SP1 , YY1 , ESR1 , p53, E2F 1 , Androgen Recep- 
tor, REL-A p65, ZNF143, NANOG, c-JUN, EGR1, ETS1 and 
NFKB p50 showed an enrichment of GO processes related to 
cell cycle and cell proliferation. We also uploaded the 
up-regulated genes from shWTl Wilms3 knockdown experi- 
ments as input list for the MetaCore program 'What are the 
key transcription factors and target genes in my data?' and iden- 
tified enriched transcription factors and related GO processes 
(Supplementary Material, Table S7). It is recognized that differ- 
ent GO processes were identified related to cell metabolism, or- 
ganismal development, differentiation and many others. 

The bioinformatics data provide good evidence that mutant 
WT1 proteins modulate transcription factor activity, and we 
wished to further explore this issue and selected p53 as a 
model system. Wild-type WT1 interacts with p53, and we 
have shown that mutant WTl Wilms2 and WTl^ 1 " 133 proteins 
also bind to p53 (Fig. 2). Wilms2 and Wilms3 cells express 
very low levels of p53 protein, which are below the detection 
limit of western blot analysis. However, we could detect serine 
15 (SI 5), serine 46 (S46) and serine 392 (S392) phosphorylation 
of p53 using an antibody array (Supplementary Material, 
Fig. S10, Panel A). To investigate whetherp53 is transcription- 
ally active, we used p5 3 -responsive and p53 non-responsive re- 
porter plasmids (pGl3-Luc and MG1 5 -Luc) that contain p53 
wild-type and p53 mutant binding sites. The relative luciferase 
activity of the p5 3 -responsive reporter compared with the p53 
non-responsive reporter was highly significant, demonstrating 
that p53 in Wilms2 and Wilms3 tumour cells is transcriptionally 
active (Supplementary Material, Fig. S10, Panel B). 

The MetaCore program 'What are the key transcription 
factors and target genes in my data?' revealed that the group of 
genes (n = 1076) that are down-regulated by shWTl contains a 
subset of 54 p53-target genes where p53 is directly engaged in 
transcriptional regulation. This information is displayed in a 
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Figure 7. Identification of p53 networks in cohorts of WT1 w,lms3 target genes. The MetaCore program 'What are the key transcription factors and target genes in my 
data' was used to identify p53 networks in subsets of genes whose expressions are either down-regulated (A) or up-regulated (B) by sh WT1 knockdown in Wilms3 cells. 
Blue circles mark down-regulated genes (A), and red circles markup-regulated genes (B). The properties of individual genes are indicated by coloured symbols. These 
symbols are explained in Supplementary Material, Table SI 5. Green arrows point to genes that are activated by p53 and red arrows to genes that are repressed by p53. 



p53 network diagram (Fig. 7A). The p53 network diagram con- 
sists of a central p53 node; red arrows are connected to p53-target 
genes whose expression is inhibited (n = 20), and green arrows 
are connected to genes whose expression is activated by p53 
(n = 34). If all genes ofthisp53 network are positively regulated 
by WT 1 Wllms3 , we can conclude that the mutant WT 1 protein can 
either enhance or antagonize the transcriptional activity of p53. 
Similar data were obtained when we investigated the group of 
genes (n = 1012) that are up-regulated by sh WT1 . Here we iden- 
tified a subset of 44 p53-target genes where 34 are activated and 
10 are repressed by p53 (Fig. 7B). These results show that the 
WT 1 Wllms3 protein has the ability to modulate the transcriptional 
activity of p53 where it most likely acts either as a co-repressor or 
co-activator in the context of individual genes. It is important to 
note that the transcriptional activity of p53 in general is controlled 
by a highly complex set of co-activators and co-repressors (30). 

To further investigate the influence of WT1 Wllms3 onp53 tran- 
scriptional activity, we used a p5 3 -responsive reporter construct 
(pGl 3-Luc) in transient transfection assays of Wilms3 cells after 
infection with shWTl lentiviruses (Supplementary Material, 
Fig. Sll). As control, a p53 non-responsive reporter construct 
(MG \5-Luc) was used in transient transfection assays. This 
initial experiment revealed diminished transcriptional activity 
of p53 in sh WT1 -treated Wilms3 cells (Supplementary Material, 
Fig. Sll, right panel). We can conclude that WT 1 Wllms3 acts as a 
co-activator for p53 in the context of this p53-responsive pro- 
moter in pG13-Zmc. 

Classification of WT1 mutations as gain-of-function 
mutations 

Next, we wished to validate the si WT1 and sh WT1 knockdown data 
from Wilms tumour cells by expressing WT 1 Wlms3 protein in a cell 
line devoid of WT1. It is a necessary prerequisite to identify the 



spectrum of genes whose expression is controlled by mutant 
WT1 in order to investigate in a second step whether the mutant 
protein gained any novel functions in comparison with the wild- 
type protein. To address these issues, we selected human MSCs 
as a model system as these do not express WT1. It is important 
for these experiments to use a cell line related to Wilms tumour 
cells in order to avoid major context-specific effects. We previous- 
ly showed that the transcriptomes and differentiation potentials of 
MSC and our Wilms tumour cells are similar (27). To express wild- 
type and mutant WT 1 in MSCs, we generated lentivirus constructs 
that express the wild-type WT1 + ~ , WT1 ++ , and the mutant 
WTl mtms3+/ - and WTr^ +l+ isoforms. After infection of 
MSCs, all four lentivirus constructs showed similar expression 
levels as determined by Q-RT-PCR (Supplementary Material, 
Fig. S12). Gene expression profiles were established using 
the same parameters as in WT1 knockdown experiments. 
In MSCs infected with wri mims3++ lentiviruses, 425 genes 
were up-regulated and 556 down-regulated. Infection with 
yyjl Wiims3+- le^v^ses ca used an up-regulation of 605 and a 
down-regulation of 755 genes. The top 100 up- and down- 
regulated genes in MSCs infected with ivTl mims3+ ~ lentiviruses 
are listed in Supplementary Material, Tables S8 and S9, respective- 
ly. We selected the subset of up-regulated genes from 
WT1 Wllms3 + " lentivirus infected MSCs and performed an enrich- 
ment analysis of pathway maps using MetaCore. The pathway 
maps 'Cell cycle_The metaphase checkpoint', 'Cell cycle_Chro- 
mosome condensation in prometaphase' and 'Cell cycle_Spindle 
assembly and chromosome separation' are most significant 

(Fig. 8A, orange bars). Similar results were obtained using 
pFn wiims3++ lentiviruses ( Fig 8A) Wue bars ) These results 

show that the mutant WT 1 protein can activate the expression of 
a significant number of genes from the G2/M phase of the cell 
cycle. A comparison of pathway enrichment data after WT1 
knockdown in Wilms3 cells (Fig. 5A) with those after 
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Figure 8. Mutant WT 1 w,lms3 controls expression of cell cycle genes from the G2/M phase in mesenchymal stem cells. (A) Enrichment of MetaCore pathway maps in 
the subset of genes whose expression is induced by WT1 m '"' s3 ++ an d WT1 wn '" s3 + " isoforms in MSC. The two isoforms of WT1 mlms3 were expressed in MSCs using 
lentiviral constructs. The length of orange bars indicates the P-value for enrichment of individual pathway maps in WT1 w,h " s3++ infected MSCs and blue bars show 
those for the wti WiI »' s +- isoform. (B) Most significant pathway map 'Cell cycle_The metaphase checkpoint': comparison of gene expression levels after sbWTl 
knockdown in Wilms3 cells and lentiviral expression of the mutant WT1 m,ms3+ isoform in MSCs. Red thermometers indicate up-regulated genes by expression 
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Figure 9. Classification of WTl nnmsz and WTr'""' J mutations as gain- 
of-function mutations. (A) Enrichment of MetaCore pathway maps for genes 
induced by wild-type WT1 in MSCs. The lengths of orange and blue bars corres- 
pond to the significance level of individual maps. Orange bars refer to WT1 ++ - 



WT1 lms overexpression in MSCs (Fig. 8A) shows complete 
concordance. This is also demonstrated by the pathway map 
'Cell cycle_The metaphase checkpoint' where the shWTl knock- 
down data from Wilms3 cells and the WT1 Wllms3+ ~ expression in 
MSCs have been compared (Fig. 8B). To illuminate quantitative 
aspects, we compared the fold change values of shWTl down- 
regulated genes (in Wilms3 cells) with those of WTl mims3+ ~- 
induced genes in MSCs, and a list of common regulated genes is 
shown in Supplementary Material, Table S 10. This gene list com- 
prises a large group of cell cycle genes from G2/M, and it is evident 
that overexpression of WT1 llms + ~ in MSCs induced these genes 
to a much higher level compared with the fold change repression by 
shWTknockdown. We conclude here that the mutant WT1 protein 
is connected to the cell cycle where it is involved in regulating ex- 
pression of genes that promote the G2/M transition and the prolif- 
eration of Wilms tumour cells. 

To address the issue whether mutant WT 1 proteins gained any 
function in comparison with wild- type WT1, we infected MSCs 
with lentiviruses expressing the wild- type WT1 ++ and WTl + ~ 
isoforms. Gene expression analysis was done as described 
earlier, and we generated datasets of up- and down-regulated 
genes. Lentiviruses expressing the WT1 ++ isoform induced the 
expression of 882 individual genes, and 690 genes showed 
reduced expression levels in MSC. Infection with WTl + ~ lenti- 
viruses was associated with an up-regulation of 425 and a down- 
regulation of 696 genes. The top 100 up- and down-regulated 
genes in MSCs infected with jYTl wihas3+ ~ lentiviruses are 
listed in Supplementary Material, Tables SI 1 and SI 2, respective- 
ly. We then performed a MetaCore pathway enrichment analysis 
using the datasets for the WT1 ++ and WTl + ~ isoforms (Fig. 9A). 
This figure shows that identical pathway maps are significantly 
enriched in gene expression profiles after infection of MSCs 
with the two WT 1 isoforms. The P- values for pathway enrichment 
are indicated for WT1 ++ and WTl + ~ by orange and blue bars, 
respectively (Fig. 9A). Although identical pathways were 
enriched, the P-values indicate subtle differences between the 
WT1 ++ and WT1 + ~ -dependent gene expression profiles. For 
the WT1 ++ isoform 'Development regulation of epithelial to mes- 
enchymal transition' and for the WT1 +I ~ isoform 'Cytoskeleton 
remodeling_Regulation of actin cytoskeleton by RHO GTPases' 
were the most significant pathway maps. Next, we looked for 
any enrichment of cell cycle pathway maps in gene expression 



induced gene expression data and blue bars to gene expression data induced by 
the WTl + ~ isoform. There was no significant enrichment of any pathway map 
related to the cell cycle. (B) Gene expression data from MSCs infected with 
]/yjl w,,ms3 + lentiviruses show proof of gain-of-function. To generate the 
heatmap, we selected a subset of genes from the cell cycle whose expression 
was down-regulated by siWTl and shWTl knockdown in Wilms tumour cells 
and compared their expression levels in diverse experiments that are indicated 
in the Figure. The colour corresponds to the log ratio of WT1 knockdown 
versus scrambled control transfected cells (blue negative log ratio and yellow 
positive log ratio). Lanes Wilms3 SH-48h. 1 and Wilms3 SH-48h.2 are biological 
replicates of lentivirus-mediated shWTI knockdown in Wilms3 cells, 48h after 
infection, and are log ratios versus scrambled RNA transduced cells. Wilms3 
SI3-48h. 1 and Wilms3 SI3-48h.2 are biological replicates of siWTl -mediated 
knockdown in Wilms3 cells 48 h after transfection versus non-silencing RNA. 
Wilms2 SI3-48h corresponds to si WT1 -mediated knockdown in Wilms2 cells. 
MSC WTl w " ms3++ and MSC WTl w " ms3+ ' refer to MSCs infected with lenti- 
viruses expressing two isoforms of mutant WTJ , and MSC wt-WTl ++ 
and MSC wt-W77 + ~ refer to MSCs infected with lentiviruses expressing two 
isoforms of wild-type WT1, expressed as log ratios versus control virus- 
transduced cells. 
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data from WT1 ++ and WT1 + -infected MSCs and visually 
inspected the top 50 enriched pathway maps (Supplementary Ma- 
terial, Fig. S13). We did not find any pathway map related to the 
cell cycle in this analysis. We can conclude from this data that 
mutant WT 1 differs from wild-type WT 1 by the capacity to regu- 
late 'cell cycle genes'. This function is novel, and therefore, these 
mutations can be classified as gain-of-function mutations. 

To confirm this result at the level of individual genes, we gen- 
erated a heat map using the down-regulated 'cell cycle genes' 
from WT1 knockdown experiments (Fig. 9B). The lanes labelled 
Wilms3 and Wilms2 show all experimental data of WT1 knock- 
down experiments using either s\WTl or shWTl lentiviruses and 
are represented as log ratios versus scrambled control RNAs. 
The lanes MSC WT1 mims3 show the log ratios of gene expression 
data after infection of MSC with lentiviruses expressing mutant 
WT1 mlmsS ) versus control virus-transduced cells, and the lanes 
MSC wt- WT1 show the gene expression levels after infection 
of MSC with lentiviruses expressing wild-type WT1 versus 
control virus. The heat map visualizes our results and shows 
that all 'cell cycle genes' that are repressed by siWTl and 
shWTl knockdown approaches are unregulated by lentivirus- 
mediated expression of mutant WTr ilms ^ in MSC. Moreover, 



lentivirus-mediated expression of wild-type WT1 in MSC 
failed to activate expression of these genes. We conclude that 
the WT1 mutations that have been identified in Wilms2 and 
Wilms3 tumour cells are gain-of-function mutations. This 
gain-of-function phenotype of mutant jypj Wllmsi [ s a j so 
evident when the top 100 genes induced by WTl Wllms3+ ~ and 
wtWTl + ~ lentiviruses are compared (Supplementary Material, 
Table SI 3). In this cohort of genes, only APOL6 and MNS1 
are induced by both lentiviruses (intensity > 200). It is very inter- 
esting to note that in the cohort of the top 100 genes that were 
down-regulated by ^77 Wilms3 +" an d wtWTl + ~ lentiviruses 
(Supplementary Material, Table S14), most genes are repressed 
by both lentivirus types. This is also visualized in the Venn dia- 
grams shown in Fig 10. Here the large overlap of genes down- 
regulated by mutant and wild-type WT1 is evident (Fig. 10A), 
whereas only a smaller set of genes is up-regulated by both len- 
tiviruses (Fig. 10B). We can conclude that the gain-of-function 
mutation in mutant Wilms tumour proteins is confined to tran- 
scriptional activation, whereas transcriptional repression is medi- 
ated by both, the wild-type and mutant WT1 proteins. Future 
work is needed to exactly work out how transcriptional activa- 
tion and repression by mutant WT1 is brought about. 



A Down-regulated genes in MSC 




B Up-regulated genes in MSC 




Figure 10. Comparison of mutant and wild-type WT1 -dependent gene expres- 
sion in MSCs. Lentiviruses expressing the + — isoforms of mutant WTl Wl,ms3 
and wild-type WT1 were used in these experiments. The expression data are 
displayed in Venn diagrams. (A) Comparison of down-regulated genes and 
(B) up-regulated genes. 



DISCUSSION 

Most WT1 mutant Wilms tumours have lost the wild-type WT1 
allele, following Knudson's two-hit hypothesis for the inactiva- 
tion of a tumour suppressor gene (31,32). In contrast in many 
adult human tumours, WT1 is highly expressed and sometimes 
mutations are found, suggesting that it acts as an oncogene 
(32,33). Some Wilms tumours with WT1 mutations express the 
mutant RNA and protein (27). Wilms tumour-derived Wilms2 
and Wilms3 cell lines were used here to address the question 
whether WT1 mutations have an impact for the biology of 
Wilms tumours. To investigate whether mutant WT1 llms2 
and WT 1 wilms3 proteins gained a new function (gain-of-function 
mutation), we compared the gene expression profiles of MSCs, 
infected with lentiviruses expressing mutant and wild-type 
WT1, and found that mutant WT1 proteins activate the expres- 
sion of a large number of genes from the G2/M phase of the 
cell cycle, whereas wild-type WT1 was unable to do so. These 
data show unequivocally that WTl Wi,ms2 and WTl Wi,ms3 muta- 
tions are gain-of-function mutations. This result is particularly 
interesting because it contradicts the currently accepted model 
that mutant WT1 proteins contribute to Wilms tumour develop- 
ment because they no longer function as tumour suppressor (32). 
Mutant WTl Wilms2 and WTl mlms3 promote the proliferation of 
Wilms tumour cells in tissue culture, and it has been shown for 
a variety of human cancers that their growth and survival can 
be impaired by the inactivation of a single oncogene, a phenom- 
enon termed oncogene addiction (34). 

This gain-of-function mutation in WT1 Wl]msi is illustrated in a 
Venn diagram (Fig. 10A and B). Here we compared the datasets 
of up-regulated genes from MSCs infected with ^; Wilms3 +- 
and wtWTl. Gain-of-function is demonstrated by the blue area 
in Figure 10B. However, the picture is more complex as there 
are genes that are regulated by both, wtWTl and the mutant 
protein (green area). In addition, as expected there are quite a 
number of genes that are induced by wtWTl only (pink area). 
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We also compared the top 100 induced genes from MSCs 
infected with wtWTl + ~ and WTl Wilms3+ ~ Antiviruses. Here it 
is evident that amongst these, only two genes were induced in 
common (Supplementary Material, Table SI 3). This clearly 
demonstrates that gain-of-function is the dominant consequence 
of this mutation. Different results are evident in the group of 
down-regulated genes (Fig. 10A). The green area indicates 
that a large fraction of genes is repressed by wtWTl + ~ and 
WT1 wttmsf + ~ i n MSCs. We also compared the top 1 00 repressed 
genes from MSCs infected with wtWTl + ~ and WTr ,,ms3+ ~ 
Antiviruses and found that the majority of genes was repressed 
by both, the wtWTl and the mutant protein (Supplementary Ma- 
terial, Table S14). Taken together these data show that a novel 
transcription activation function is the major phenotype of the 
WT1 mutation in Wilms3 cells. The repression function of 
wtWT 1 and the mutant protein may be explained by the presence 
of the Pro/Gln-rich transregulatory and dimerization domain in 
the N-terminal domain. This domain also contains a putative 
RNA recognition motif (5,8,17), and mutant Wilms3 protein 
has retained its ability to bind to RNA aptamers, albeit with 
reduced affinity. Therefore, the mutant proteins can interact 
and modulate the function of other proteins and/or bind to 
RNA as well. However, the mechanisms by which the mutant 
proteins exert these functions remain to be identified. 

It is interesting to note that a leukaemia-associated mutant WT 1 
protein devoid of the ZF domain, WT 1 (delZ) , has also gained new 
functions (35). Whereas WT 1 inhibited the proliferation of haem- 
atopoietic progenitor cells (36), the WT l(delZ) protein promoted 
expansion of human haematopoietic progenitor cells in tissue 
culture. WT1 mutations were detected in 10-15% of adult AML 
cases at diagnosis, and the presence of WT1 mutations is an inde- 
pendent predictor of poor clinical outcome (35). Gene expression 
profiling revealed that WTl(delZ) expression was associated 
with increased gene expression levels of certain genes related to 
proliferation and decreased expression levels of genes related to 
myeloid differentiation. 

Mutant WT1 controls genes from the G2/M phase of the cell 
cycle. For example, the expression of BIRC5, BUB1, CCNB1, 
CCNB2, CDC20, CDK1, CENPF, PLK1 and UBE2C was 
down-regulated by shWTl in Wilms3 cells and up-regulated in 
MSCs infected with Antiviruses expressing the WTl Wilms3+ ~ 
isoform. It is remarkable that CCNB1, CCNB2 and CDK1 are 
amongst these genes as they are key regulators of cell cycle pro- 
gression in G2/M phase. The most significant transcriptional sig- 
nature of cancer cells and malignant tumours is the up-regulation 
of cell cycle and proliferation-related genes (37). An analysis of 
gene expression data from a large variety of human cancer types 
identified a core proliferation cluster (core PC) predominantly 
consisting of genes from the G2/M phase of the cell cycle (38). 
A comparison of our gene expression data with those from the 
core PC showed a high degree of overlap. The core PC consists 
of 46 genes and 2 1 of these genes showed reduced expression in 
lentiviral shWTl knockdown experiments. Similar results were 
obtained by s\WTl knockdown in Wilms2 and Wilms3 cells. 
We can conclude from these data that the WT1 gain-of-function 
mutations in Wilms2 and Wilms3 cells are oncogenic events 
affecting gene expression of important cell cycle regulators and 
in consequence increased proliferation of Wilms tumour cells. 
This is supported by our WT1 knockdown experiments, which 
showed that the mutant WT1 proteins promote proliferation of 



Wilms2 and Wilms3 tumour cells. This finding provides convin- 
cing experimental evidence that the WT1 mutations may contrib- 
ute to the development of Wilms tumours in vivo. 

A large number of cell cycle regulatory genes are also princi- 
pal targets of p53-mediated repression (38,39). Transcriptional 
repression mediated by p53 is complex and involves other tran- 
scription factors like E2F1 and NFY (38). Mutant WT1 proteins 
can activate the transcription of cell cycle genes that are 
repressed by p53. These include CCNB1, CDK1, PLK1, 
CDC20, MAD2a and others (39-43). The promoter of the 
Cyclin B (CCNB1) gene for example contains a p53 head-to-tail 
response element, and binding of a p53 dimer to this response 
element is sufficient to abrogate transcription (40). We show 
here that mutant WT1 can activate CCNB1 and CCNB2 expres- 
sion, and this suggests that the mutant protein can antagonize the 
repression function of p53 at the Cyclin B promoter. Future 
experiments will elucidate how the activation of the Cyclin Bl 
gene expression is brought about. WT1 interacts with p53 and 
modulates its function (22). Mutant WTl wilms2 and WTl wilms3 
proteins can bind to p53 in vitro, and it is possible that this 
property enables WTl Wllms3 to neutralize p53-dependent 
repression of the Cyclin B promoter in vivo. The major pheno- 
type of the WT1 mutation is the gain of a transactivation func- 
tion, and an interaction of mutant WT1 with p53 at the Cyclin 
Bl promoter might enable this transactivation function to over- 
come repression by p53. Another possible scenario is the seques- 
tration of p53 by mutant WT1 in the cytoplasm, resulting in 
loss of p53 function in the nucleus. In the context of this issue, 
it is important to study in future experiments, whether mutant 
WT1 proteins can be detected in the nucleus, even in minor 
amounts. Our immunofluorescence analysis (Supplementary 
Material, Fig. SI) and recent work (27) showed a predominant 
cytoplasmic localization of mutant WT1 proteins. Despite of 
this, it is possible that small amounts of these proteins are 
transported to the nucleus by a piggyback mechanism after asso- 
ciation with other nuclear proteins. This is a speculative argu- 
ment at this stage; however, a mutant WT 1 protein lacking the 
ZF domain (WTl-delZ) was identified in nuclei of transfected 
osteosarcoma cells (44). 

We have shown that Wilms2 and Wilms3 proteins are able to 
interact with RNA. It is possible that cytoplasmic localization of 
mutant WT 1 proteins affects gene expression by altering the me- 
tabolism of RNA in the cells. This is, however, a complex issue 
and investigations addressing the functions of cytoplasmic 
mutant WT1 proteins are beyond the scope of the current 
work. In addition, we provided evidence based on an analysis 
using the MetaCore program 'What are the key transcription 
factors and target genes in my data? ' that mutant WT 1 modulates 
the activities of selected transcription factors also in the context 
of cell cycle-related gene expression. Again these results 
revealed an additional level of complexity, and this will be 
addressed in future studies. 

Our WT1 knockdown data show that mutant WT 1 inhibits the 
expression of a large number of genes from signalling pathways 
including several members of the canonical Wnt and planar cell 
polarity signalling pathways. The regulation of these pathway 
activities is complex, and in some cases, gene expression 
levels of inhibitors (e.g. AXIN2) can be up-regulated in the 
active state of the pathway. To ultimately determine activity of 
a given pathway of interest, functional studies are required, 
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which are beyond the scope of this investigation and will be part 
of future work. With regard to Wnt signalling, it is important to 
note that wild-type WT 1 inhibits Wnt signalling because it com- 
petes with TCF4 for binding to CBP (45). Moreover, depletion of 
wild-type WT1 results in up-regulation of CTNNB1 expression 
and WT1 negatively regulates Wnt signalling in Sertoli cells 
(46). Mutant WT1 Wilmsf down-regulates CTNNB1 expression 
that might also inhibit the Wnt signalling pathway. In this 
context, it is of special interest to note that Wilms3 cells 
express wild-type CTNNB1. Therefore, Wilms3 tumour cells 
can survive without CTNNB1 mutations that are a common 
phenotype of Wilms tumours where WT1 is mutated. Mutant 
WTl Wllms3 was classified as a gain-of-function mutation with 
oncogenic potential, and the WT1 knockdown data showed 
that the growth and survival of Wilms3 tumour cells is dependent 
on the expression of WT1 wilms3 . 

Furthermore, several genes that are repressed by mutant WT 1 
are involved in chromatin modification and some are positive 
and negative regulators of MAPK signalling. It was reported that 
wild-type WT1 negatively regulates the RAS/MAPK signalling 
pathway by inducing the expression of MKP3 (DUSP6), a 
phosphatase and negative regulator of this pathway (47). Here 
we found that mutant WT1 inhibits the expression of the dual spe- 
cificity phosphatase genes DUSP3 and D USP1 6. Another group of 
genes that are down-regulated by mutant WT1 are differentiation 
regulating genes, e.g. GREM2, HMOX1, MY05a, MY09B, 
GDNF, NES, NKX3-1, PODXL, HGFL and SNPO. By using 
mouse models, it was shown that PODXL is a WT1 target gene 
with high expression levels in podocytes and highly differentiated 
kidney epithelial cells (48). PODXL gene expression is induced by 
WT1 and repressed by p53 (49). Wild-type WT1 activates PODXL 
by direct binding to its promoter, and WT1 and PODXL are 
co-expressed in the developing kidney (50). We have reported 
recently that WT1 and BASP1 interact with the promoter of 
PODXL, BASP1 inhibiting WT1. During differentiation, BASP1 
promoter binding is down-regulated and in consequence WT1 
can activate PODXL1 gene expression (51). Since the expression 
of PODXL is inhibited by WTl mlms3 and PODXL expression is 
increased during kidney differentiation, one could argue that 
mutant WT 1 keeps the tumour cells in an undifferentiated or less 
differentiated state. 

In summary, we have identified novel gain-of-function muta- 
tions of WT1 that occurred in vivo during the development of 
Wilms tumours in patients. These mutations may drive tumour 
development as they promote cell proliferation in Wilms 
tumour cells. This result is of special interest for the biology 
of Wilms tumours as an ablation of W77 alone does not induce 
Wilms tumours in transgenic mice. However, in a transgenic 
mouse model with carrying a WT1 truncation mutation (396 
stop), one animal developed a WT (52). It is interesting that in 
this Wilms tumour, one allele carried the targeted mutation 
and the other allele showed an exon 9 skipping event. Therefore, 
this tumour developed when both WT1 alleles carried mutations, 
corresponding to our Wilms tumour cell lines, where both alleles 
harbour the same the mutation (52). These data are reminiscent 
of the situation with another tumour suppressor, p53. A mouse 
model for pancreatic cancer revealed that a mutant form of p53 
had an increased tumorigenicity through a gain-of-function 
mechanism. In this case, mutant p53 had extra functions over 
and above the loss of p53 in a tumour model (53). 



MATERIALS AND METHODS 

Cell culture 

Human mesenchymal stem cells were purchased from Lonza and 
maintained in MSCGM (Lonza). The primary Wilms tumour 
cell lines Wilms2 and Wilms3 were described previously (27) 
and were cultivated in MSCGM with a medium change every 
3 days. 

Plasmids and lentiviral expression vectors 

Full-length human wild-type and mutant (p.V432fsX87) WT1 
genes with an optimal Kozak sequence in front of the start ATG 
were cloned through PCR from mutant WTl m,ms3++ and 
WT1 mlm3+ ~ isoforms inpCRII-TOPO vectors and from wild-type 
WT1 ++ and WTl + ~ isoforms in pDNA3.3-Topo-TA vectors. 
Primers for the PCR were Kozak- ^77-for (5' AACGAATTCAC 
CACCATGGGCTCCGACGTGCGG 3') and wt- WTl-rev (5' GT 
TGAATTCTCAAAGCGCCAGCTGGAGTT 3') or mut-^77- 
rev (5' CAGTCGACTCAACTAAAAGTAGGCAGGGC 3')- 
The PCR products were cloned into the pCRII-TOPO vector and 
excised by digestion with EcoRI (wild-type W 77) or KpnI/Xbal 
(mutant WT1) and ligated into the pENTR4 vector [Addgene, 
17 424 (54)], digested with the same restriction enzymes. The 
wild-type and mutant WT1 isoforms were finally cloned into the 
pLenti CMV Puro Dest vector [Addgene, 17 452 (54)] through 
the Gateway Technology as described in the 'Gateway LR 
Clonase II Enzyme Mix' protocol (Invitrogen). As a control for 
the transduction of expression vectors, a FOP luciferase pLenti 
XI dest puro vector was used, containing the FOP expression 
cassette, cloned from the M51 super 8 FOP Flash [Addgene 
(55)]. All vector constructs were verified by sequencing. 

The construction of a plasmid vector used to express the ZF 
domain of wild-type WT1 has been described previously (56). 
The cDNA used in this study to express the ZF domain ( — KTS) 
of Wilms3 was constructed by the polymerase chain reaction 
(PCR) using an upstream primer at the start of ZF 1, containing 
a Ncol recognition sequence and a downstream primer including 
the termination codon for the Wilms3 reading frame, containing a 
EcoRI recognition sequence. PCR products were digested with 
Ncol and EcoRI and ligated into plasmid pET30a, digested with 
the same restriction enzymes. The Ncol site puts the Wilms3 
cDNA in frame with the vector-included his-tag. Plasmids were 
isolated from appropriate colonies, and the correct plasmid con- 
structs were confirmed by DNA sequencing. 

The 'Mission shRNA Plasmid DNA' with a pLKO.l-puro 
backbone used in the WT1 knockdown experiments and the 
non-mammalian shRNA control used in all transduction 
experiments were purchased from Sigma-Aldrich; shWTl 
(exon7, TRCN0000040067) and sh scrambled (SHC002). 

Production of lentiviral viral particles, transduction 
and siRNA transfection 

For virus production, HEK293 T cells were transfected with 6 |xg 
of each plasmid DNA: packaging vector pczVSV-G (57), 
pCD NL-BH (57) and the desired expression vector in the pres- 
ence of 45 |xg PEI (Polyethyleneimine, branched, Aldrich). As 
internal transfection control, a lentiviral GFP-expression Vector 
(pCL7EGwo, obtained from Dr H. Hanenberg) was transfected 
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in each experiment. The medium was changed after 1 day to 
IMDM with 10% FCS and 1% penicillin/streptomycin, and 48 h 
after transfection, viral supernatants were harvested, filtered 
(0.45 Jim), concentrated by centrifugation and cryo-conserved. 

WT1 knockdown experiments were performed with (1) trans- 
fection with siRNA oligonucleotides {siWTl) and (2) infection 
with lentiviruses expressing shRNA (sh WT1). For si WT1 knock- 
down, Wilms2 and Wilms3 cells were seeded in 24-well plates 
(4 x 10 4 cells/well), transfected after 24 h with the 'HiPerfect 
transfection reagent' (Qiagen) with a combination of two differ- 
ent siRNAs (Qiagen: Hs_WTl_l_HP in 3' UTR and 
Hs_WTl_8_HP in exon 6) and a non-silencing control siRNA, 
Alexa Fluor 488 labelled (Qiagen). Total RNA was isolated at 
the indicated time points. 

For shWTl experiments, Wilms3 cells were seeded in a 6-well 
plate (5 x 10 4 cells/well). After 24 h, the cells were infected 
with shWTl lentiviruses or as control with lentiviruses expres- 
sing a non-silencing shRNA in the presence of polybrene 
(5 |JLg/ml, H9268, Sigma), and the medium was changed after 
4 h. RNA was isolated 48 h after transduction. All experiments 
were performed in duplicates. 

Overexpression of wild-type and mutant proteins in MSCs 
was achieved by seeding 1 x 10 5 MSC in 6-well plates and in- 
fection after 24 h with 2 ml of the respective lentiviruses in the 
presence of polybrene (1 mg/ml). For a biological replicate, 
the experiments were done in duplicates. The medium was 
changed after 24 h to MSCGM, and RNA was isolated at 48 h. 
The negative control was a scrambled lentiviral vector. 

RNA isolation, gene expression analysis and Q-RT-PCR 

Total RNA was isolated after WT1 knockdown or overexpres- 
sion using the RNeasy Mini Kit (Qiagen). Labelling of total 
RNA in the Two-Colour format was performed as described 
by the manufacturer (Agilent Technologies). The cRNA was 
hybridized to 4 x 44 K 'Whole Human Genome Oligo Microar- 
rays' (Agilent, VI), and 'Spike-In positive controls' (Agilent) 
were added as internal controls. Expression arrays were 
scanned (Agilent Technologies Scanner G2505B) and quantified 
using the Agilent Feature Extraction V10. 1.1.1 Software. The 
software package LIMMA (R-library, bioconductor, http 
://www.R-project.org) was used for statistical analysis and the 
search for differentially expressed genes. The gene expression 
data were submitted and are found at GEO, GSE54635. 

For Q-RT-PCR analyses, cDNAs were synthesized using 
TaqMan Reverse Transcription Reagents (Applied Biosystems). 
All Q-RT-PCR experiments were performed in triplicates using 
the TaqMan gene expression assay probes from Applied Biosys- 
tems for WT1 (Hs00240913_ml),M£'SDC2(Hs00398059_ml), 
GDF5 (Hs00167060_ml) CASP7 (Hs01029847-_ml), PODXL 
(hCG18687) and JARID2 (Hsl004460_ml) with the TaqMan 
Universal PCR master Mix (no AmpEraseUNG) or Brilliant II 
QPCR Master Mix with Low Rox (Agilent Technologies). In 
the siWTl experiments, the expression levels were normalized 
with the 7&SRNA (Hs99999901_sl), in the shWTl experiments 
with GAPDH (4333764 T). The WT1 knockdown was normal- 
ized versus ACTB (4333762F) and the overexpression versus 
RER1 (Hs00199824_ml). The reason for using these different 
calibrator probes was that we noted a regulation of the calibrator 
gene in some of these experiments, and therefore, we have 



identified a new gene that was basically never regulated in any 
of our gene expression studies: RER1. The siWTl Q-RT-PCRs 
were run on an ABI PRISM 7900HT, and the Q-RT-PCRs 
from shWTl and expression of WT1 in MSC were run on a 
Mx3000P Sequence Detection System (Stratagene). 



Microarray data analysis 

Microarray data analysis was performed using the MetaCore 
software suite (http://thomsonreuters.com/metacore/) to iden- 
tify biological processes associated with gene expression pro- 
files. For analysis of gene expression profiles, normalized 
microarray data were manually curated before uploading indi- 
vidual data files. We selected only named genes from datasets 
and used a cut-off of 200 fluorescence intensity and a fold 
change of 1.5. The Agilent Identification numbers of selected 
genes were then uploaded and analysed using the MetaCore soft- 
ware program Enrichment Analysis (EA) for Datasets: 

Enrichment analysis consists of mapping gene IDs of the data- 
sets) of interest onto gene IDs in entities (terms) of built-in func- 
tional ontologies represented in MetaCore™/MetaDrug™ by 
pathway maps and networks. Mapping procedure involves cal- 
culating statistical relevance of the matches found. We use 
both ontologies, which are the result of our in-house extensive 
classification effort, and public ontologies such as Gene Ontol- 
ogy (GO, www.geneontology.org). Ontologies are represented 
by canonical pathway maps, cellular process networks, disease 
biomarker networks, drug target networks, toxicity networks 
and metabolic networks. GO ontologies include GO processes, 
GO molecular functions and GO localizations. 

Statistical relevance of the found ontology matches is calcu- 
lated as P-value, or a probability of a match to occur by 
chance, given the size of the database. The lower the P-value, 
the higher is the 'non-randomness' of finding the intersection 
between the dataset and the particular ontology. That, in turn, 
translates into a higher rating for the term matched. In particular, 
the more genes/proteins belong to a process/pathway, the lower 
the P-value. MetaCore program 'What are the key transcription 
factors and target genes in my data? ' uses a database containing 
all published information regarding transcription factors and 
their target genes. We uploaded the up- or down-regulated 
genes from shWTl knockdown experiments in Wilms3 cells as 
the input list. The transcription regulation workflow generates 
sub-networks centred on transcription factors. Sub-networks 
are ranked by a P-value and interpreted in terms of Gene Ontol- 
ogy. The resulted sub-networks of the whole network are dis- 
played in Supplementary Material, Table S6 and Table S7, a 
network list table, in the order reflecting their statistical signifi- 
cance in terms of P-values. 



Cell proliferation after WT1 knockdown 

For the determination of the proliferation rates of Wilms2 and 
Wilms3 cells after WT1 knockdown, 5 x 10 4 cells per well 
were seeded in MSCG medium in a 6-well plate. After transduc- 
tion with sliWTl and scrambled control lentiviral particles, the 
cells were counted in a Neubauer chamber 24 h, 48 h, 72 and 
96 h. Each time point was analysed in duplicate. RNA was 
extracted from these cells for gene expression analysis. 
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Expression and purification of recombinant ZF proteins 

Preparations of N-terminal his-tagged ZF proteins were carried 
out as described previously (58). Protein purity was confirmed 
by SDS-polyacrylamide gel electrophoresis (PAGE), and the 
concentration of each protein preparation was determined by 
the method of Bradford. 



Radiolabelling of RNA ligands and equilibrium binding 
of RNAs to ZF proteins 

Preparation of radiolabelled RNA aptamers by in vitro transcrip- 
tion was as described previously (29). The apparent association 
constants for the binding of radiolabelled RNAs to ZF proteins 
were determined using a double-filter binding assay. The 
binding buffer consisted of 20 mM Tris-HCl pH 7.5 (20°C), 
5 mM MgCE;, 100 mM KC1, 10 |xm ZnCL?, 0.5 mM Tris (2-carbox- 
yethyl)phosphine hydrochloride, 100 |Jig/ml BSA and 1 |xg/ml 
polydldC. The affinity of each ZF protein for RNA was determined 
using three or more independent assays. Apparent dissociation 
constants (K d ) for the binding of the mutant and wild-type proteins 
to RNA were calculated by fitting the data to a simple bimolecular 
equilibrium model using the general curve fitting function of 
Kaleidagraph software (Synergy Software, Reading, PA, USA) 
and the following equation: 

[RNA-protein] _ [protein] total 
[RNAJtotai [protein] total + K A 

where [RNA] total << K d and [RNA-protein]/[RNA] total is 
reported as the fraction of RNA bound. Association constant 
(K a ) values were derived as the reciprocal of the measured K d 
values and are reported as the mean of three independent determi- 
nations with the associated standard deviations. Relative affinities 
were arrived at by dividing the apparent K a for the mutant protein 
by the apparent K a for the wild-type protein determined in parallel, 
or by dividing the apparent K & for the mutant RNA by the apparent 
K & for the wild-type RNA determined in parallel. The errors for 
relative affinities are given by the expression cr= {(cr 1/M1) 2 + 
(cr2/M2) 2 } 1/2 x M2/M1, where Ml and M2 are the respective 
association constants for wild-type and mutant protein and the 
cr values are the corresponding standard deviations for these 
determinations. 



Analysis of endogenous p53 activity in Wilms cells 

Endogenous p53 activity in Wilms3 cells was measured with 
the p53-responsive reporter pG13-Z,»c [13 p53 binding sites, 
Addgene, 16 442 (59)] and the non-responsive reporter 
pG!5-Luc [15 mutant p53 binding sites, Addgene, 16443 (59)]. 
As control, Wilms3 cells were transfected with a promoter-less 
pGL3-basic luciferase vector. All cells were co-transfected with 
pGL4.75 (Renilla-Luc) for normalization. Cells were lysed after 
24 h in 1 x Passive Lysis Buffer (Promega), and the luciferase 
activity was measured with the Beetle-/Renilla-Juice Big Kit 
(PJK GmbH, Germany) as described by the manufacturer. The 
relative p53 luciferase activity was determined as quotient of 
the firefly-luciferase and renilla-luciferase activity. The analysis 
was performed in triplicates. 



For the determination of the status of the p53 protein, cell lysates 
from Wilms2 and Wilms3 cell lines were used for probing the 
Human Phospho-Kinase Array (Catalog # ARY003B, R&D). 
The analysis was performed as described by the manufacturer 
(R&D). 



Protein interaction of WT1 and p53 

HeLa cells were cultured in Dulbecco ' s modified Eagle ' s medium 
and transfected with Effectene (Qiagen) with pAcGFP vectors 
expressing N-terminally GFP-tagged WTl wilms3 or WTl Witals2 
or wild- type WT1 fusion proteins (27). Preparation of whole 
cell extracts and IP western blot were carried out as described pre- 
viously (60). Antibodies anti-WTl (F6) and anti-p53 (DOl) 
sc-126 were obtained from Santa Cruz Biotechnology. 
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